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We  present  a  magnetophotoluminescence  study  of  individual  vertically  stacked  InAs/GaAs  quantum 
dot  pairs  separated  by  thin  tunnel  barriers.  As  an  applied  electric  field  tunes  the  relative  energies  of  the  two 
dots,  we  observe  a  strong  resonant  increase  or  decrease  in  the  g  factors  of  different  spin  states  that  have 
molecular  wave  functions  distributed  over  both  quantum  dots.  We  propose  a  phenomenological  model  for 
the  change  in  g  factor  based  on  resonant  changes  in  the  amplitude  of  the  wave  function  in  the  barrier  due 
to  the  formation  of  bonding  and  antibonding  orbitals. 
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Quantum  dots  and  quantum  dot  molecules  (QDMs)  have 
proven  to  be  a  versatile  medium  for  isolating  and  manipu¬ 
lating  spins  [1,2],  which  are  of  great  interest  for  quantum 
information  processing  [3,4].  In  particular,  photolumines¬ 
cence  (PL)  spectra  have  been  used  in  self-assembled 
QDMs  to  observe  coherent  tunneling  [5-8]  and  identify 
spin  interactions  through  fine  structure  [9].  Electrical  con¬ 
trol  of  isolated  spins  through  their  g  factors  is  highly 
desirable  for  implementation  of  quantum  gate  operations. 
To  date,  electrical  control  of  g  factors  has  only  been 
observed  in  ensembles  of  electrons  in  quantum  wells  by 
shifting  the  electron  wave  functions  into  different  materials 
[10-13].  In  this  Letter  we  present  a  striking  electric  field 
resonance  in  the  g  factor  for  molecular  spin  states  confined 
to  a  single  quantum  dot  molecule. 

To  our  knowledge  this  is  the  first  observation  of  electri¬ 
cal  control  over  the  g  factor  for  a  single  confined  spin. 
Moreover,  the  isolation  of  a  single  QDM  allows  us  to 
spectrally  resolve  and  identify  individual  molecular  spin 
states  that  have  different  g  factor  behaviors.  In  Fig.  1(a)  we 
indicate  molecular  spin  states  of  both  the  neutral  exciton 
(Z°,  one  electron  recombining  with  one  hole)  and  positive 
trion  (Z+,  electron-hole  recombination  in  the  presence  of 
an  extra  hole)  at  zero  magnetic  field.  The  different  electric 
field  dependences  of  the  g  factors  for  these  states  is  appar¬ 
ent  in  Fig.  1(b),  where  the  splitting  of  PL  lines  increases  for 
some  molecular  spin  states  and  decreases  for  others.  This 
electric  field  dependence  is  nearly  an  order  of  magnitude 
larger  than  previously  reported  in  quantum  wells  [10-13]. 
The  effect  arises  from  the  formation  of  bonding  and  anti¬ 
bonding  orbitals,  which  results  in  a  change  in  the  ampli¬ 
tude  of  the  wave  function  in  the  barrier  at  resonance. 

Our  QDMs  consist  of  two  vertically  stacked  self- 
assembled  InAs  dots  truncated  at  a  thickness  of  2.5  nm 
and  separated  by  a  2  nm  GaAs  tunneling  barrier  [  14].  As  an 
applied  electric  field  tunes  the  relative  energies  of  the  two 


dots,  strong  tunnel  coupling  between  the  hole  states  creates 
the  molecular  spin  states.  Unlike  samples  with  a  thicker 
tunnel  barrier  [5],  the  states  retain  molecular  character 
throughout  the  observed  range  of  electric  fields.  We 
present  data  from  a  single  molecule,  but  the  universality 
of  the  behavior  has  been  verified  by  detailed  studies  of 
7  other  molecules  from  the  same  sample.  We  first  explain 
the  spectra  and  molecular  spin  states  at  B  =  0  T.  We  then 
describe  the  magnetic  field  dependence  and  propose  a 
phenomenological  model  for  the  electric  field-dependent 
Zeeman  splitting. 

In  Fig.  2  we  show  all  PL  lines  from  X°  and  X+  at  B  = 
0  T.  These  lines  are  identified  by  their  relative  energies,  the 
power  dependence  of  their  intensities  and  the  electric  field 
dependence  of  the  anticrossings  [5],  The  X{)  lines 
[Fig.  2(a)]  show  a  clear  anticrossing  at  Fx The  anticross¬ 
ing  arises  from  tunnel  coupling  between  the  ground  state 
hole  levels  in  each  dot,  which  forms  bonding  (X[]ls)  and 
antibonding  (Y({)  orbitals.  The  electron  remains  in  the 
bottom  dot  throughout  the  range  of  electric  fields  consid¬ 
ered  here  [5].  The  X+  lines  [Fig.  2(b)]  have  a  more  com- 
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FIG.  1  (color  online),  (a)  B  =  0  T  photoluminescence  spectra 
of  a  single  QDM.  The  complex  pattern  of  anticrossings  arises 
from  the  formation  of  molecular  spin  states,  (b)  At  B  =  6  T,  the 
molecular  spin  states  have  a  Zeeman  splitting  (bars)  that  depends 
strongly  on  the  applied  electric  field  (F). 
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FIG.  2  (color  online).  Energies  extracted  from  Fig.  1(a).  (a)  X° 
lines  anticross  at  Fx o ,  where  the  direct  (lower  inset)  and  indirect 
(upper  inset)  transitions  are  degenerate,  (b)  X+  initial  states 
anticross  at  Fx+  and  near  90  kV/cm.  Final  (  hole)  states  anticross 
at  Fh.  Inset:  charge  distribution  for  the  circled  singlet  and  triplet 
transitions. 


plicated  pattern  because  anticrossings  occur  in  both  the 
initial  (one  electron  and  two  holes)  and  final  (one  hole) 
states  [5,9]. 

To  explain  the  origin  of  the  X+  molecular  spin  states  we 
turn  to  the  Hamiltonians.  The  basis  states  will  be  identified 
as  where  eB  [eT]  are  the  spins  of  electrons  (±1/2: 

t  or  | )  and  hB  [hT]  the  spins  of  heavy  holes  (±3/2:  K  or  jj. ) 
in  the  bottom  [top]  dot.  X  indicates  an  exciton  (h  a  single 
hole)  and  Q  is  the  net  charge,  k  is  the  total  spin  projection. 
Singlets,  which  have  total  spin  ±1/2,  will  be  denoted  Xs  to 
distinguish  them  from  the  ±1/2  triplets. 

The  final  state  has  only  a  single  hole  with  relative 
energies  given  by  the  diagonalization  of  the  Hamiltonian: 


Hh  = 


0  h  \ 

th  edF  j 


(1) 


in  the  spin-degenerate  {jo^+3/2>  o’tr^+3/2  basis,  where  th  is 
the  tunneling  energy  for  a  bare  hole,  d  is  the  distance 
between  dot  centers  and  F  is  the  applied  electric  field. 
The  energies  of  the  final  states  as  a  function  of  electric 
field  are  shown  in  Fig.  3(b).  The  formation  of  bonding  and 
antibonding  states  at  the  anticrossing  point  is  illustrated. 

The  initial  state  (X+)  contains  one  electron  and  two 
holes.  For  simplicity  we  present  only  the  electron- spin-up 
case,  which  is  degenerate  with  the  spin-down  case  at  zero 
magnetic  field.  If  both  holes  are  in  the  same  dot,  the  Pauli 
principle  requires  singlet  configurations:  }||/qX/  (both  in 
the  bottom  dot)  or  (both  in  the  top  dot).  If  there  is 

one  hole  in  each  dot,  singlet  (J^X/-)  and  triplet  (j'jjw , /2, 
jj°X/7/2)  configurations  are  possible.  By  J^X^  we 
mean  the  antisymmetric  hole  spin  wave  function  ( jj. ,  1) 
-tr.iuiX  , j2  means  the  symmetric  hole  spin  wave 
function  ( j) ,  ft  +  ]]',]].).  The  relative  energies  of  the  initial 
states  are  given  by  the  Hamiltonian 


Electric  Field  (kV/cm) 


FIG.  3  (color  online).  Calculated  zero  magnetic  field  energies 
of  spin-degenerate  (a)  initial  states  (X+)  and  (b)  final  states  (h+) 
[24] .  Labels  indicate  the  dominant  basis  state  for  easy  compari¬ 
son  to  Eq.  (2).  Insets  in  (b)  show  the  bonding  (lower)  and 
antibonding  (upper)  wave  functions  (green)  and  probability 
distributions  (red),  (c)  Fine  structure  of  singlet  and  triplet  states 
at  F  =  79.2  kV/cm.  Arrows  indicate  optically  allowed  transi¬ 
tions  to  the  bonding  final  state. 
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basis.  I/  and  T2  are  due  to  Coulomb  interactions,  tx+  is 
the  Coulomb-corrected  tunneling  energy  for  a  hole  in  the 
presence  of  an  electron  and  additional  hole,  Ek  =  edF  ± 
mJeh  with  m  =  0,  +1,  -1  for  k  =  ±1/2,  -5/2,  ±7/2. 
Jeh  is  the  exchange  energy  between  an  electron  and  hole  in 
the  same  dot  [5,15],  Diagonalizing  Eq.  (2)  gives  the  en¬ 
ergies  of  the  initial  states,  which  are  plotted  as  functions  of 
electric  field  in  Fig.  3(a). 

Because  tunneling  is  a  spin  conserving  process,  only  the 
j^X/  singlet  state  can  tunnel  couple  with  l/^X/  and 
})7:jX/,  which  must  be  singlets  because  the  two  holes  are 
in  the  same  dot.  These  3  singlets  are  therefore  strongly 
mixed  to  create  molecular  orbital  states  that  anticross  near 
57  kV/cm  (Fx+)  and  90  kV/cm.  Unlike  the  singlet  states, 
the  triplet  states  do  not  mix  and  are  not  affected  by  these 
anticrossings  [16].  This  creates  a  “kinetic”  splitting  be¬ 
tween  triplet  and  singlet  states  [17].  An  example  is  indi¬ 
cated  by  the  dashed  oval  in  Fig.  2(b),  where  the  lower 
energy  singlet  line  remains  separated  from  the  two  opti- 
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cally  allowed  triplet  lines,  which  are  split  by  electron-hole 
exchange.  The  fine  structure  of  the  corresponding  states  is 
shown  in  Fig.  3(c). 

By  adding  the  Zeeman  interaction  to  the  Hamiltonian 
[18],  we  calculate  the  magnetic  field  dependence  of  the 
molecular  spin  states.  Figure  4(a)  shows  the  states  of 
Fig.  3(c),  which  split  into  doublets  with  an  applied  longi¬ 
tudinal  magnetic  field.  The  final  states  are  simply  the  two 
spin  orientations  of  a  single  hole,  split  by  gh.  The  splitting 
of  initial  states  depends  on  their  spin  configuration. 
Because  of  the  two  parallel  hole  spins,  has  a  large 

splitting  given  by  ge  +  2 gh  while  X+7/2  is  split  by  -ge  + 
2 gh.  In  contrast,  the  ±1/2  singlet  (V/")  and  triplet  ( WT ,  /2) 
have  oppositely  paired  hole  spins  and  therefore  a  small 
splitting  given  by  ge.  The  g  factor  for  PL  transitions  is 
given  by  the  difference  in  g  factor  between  the  initial  and 
final  states.  Away  from  the  electric  field  resonances,  the  g 
factor  is  gT  =  ge  +  gh,  as  indicated  by  the  vertical  lines. 

Using  the  model  described  below,  we  obtain  gT  =  —2.2. 
To  plot  the  initial  and  final  states  in  Fig.  4(a)  we  have  taken 
the  relative  weights  of  the  electron  and  hole  g  factors  to 
match  those  obtained  by  Bayer:  ge  =  —0.6  and  gh  =  —  1.6 
[18].  The  calculated  energies  of  the  transitions  are  shown 
by  the  lines  in  Fig.  4(b).  The  experimentally  observed  PL 
energies  are  given  by  the  symbols,  with  the  diamagnetic 
shift  (10.9  /reV/T2)  subtracted. 

The  g-factor  resonances  are  clearly  evident  in  Fig.  5, 
where  the  symbols  plot  the  measured  energy  splitting  of 
the  X°  and  two  X+  Zeeman  doublets  at  B  =  6  T  as  a 
function  of  electric  field.  Strong  enhancement  or  suppres¬ 
sion  of  the  splitting  is  observed  at  Fh  and  Fx o.  All  of  the 


FIG.  4  (color  online),  (a)  Magnetic  field  dependence  of  initial 
(X+)  and  final  (/i)  states  for  the  singlet  and  triplet  transitions  at 
an  electric  field  of  79.2  kV/cm  [schematic  Fig.  3(c)].  Vertical 
lines  indicate  spin  allowed  recombinations,  (b)  Calculated  (col¬ 
ored  lines)  and  experimentally  observed  (black  points)  PL  en¬ 
ergies.  (c)  Zeeman  splitting  of  PL  lines  (A)  and  initial  and  final 
states  for  Xt$/2  (red  lines)  and  (green  lines). 


data  can  be  qualitatively  explained  by  a  phenomenological 
model  of  the  formation  of  bonding  and  antibonding  orbi¬ 
tals,  which  results  in  resonant  changes  in  the  amplitude  of 
the  wave  function  in  the  barrier.  We  focus  first  on  transi¬ 
tions  involving  a  bonding  orbital  (open  symbols),  which 
have  an  increased  amplitude  of  the  wave  function  in  the 
barrier  [Fig.  3(b)  lower  inset]. 

The  wave  function  for  the  bonding  orbital  of  a  single 
hole  can  be  written  as  =  a\  1)  +  b\2),  where  1 1)  and  |2) 
are  the  basis  states  of  Eq.  (1),  the  wave  functions  for  holes 
localized  in  the  two  different  dots.  The  coefficients  a  and  b 
are  functions  of  electric  field  determined  by  Eq.  (1).  The 
electric  field-dependent  g  factor  for  a  hole  in  a  bonding 
orbital  is  given  by  gf  (F)  =  (^Ig^z)!^},  where  gh(z)  is 
the  hole  g  factor  as  a  function  of  position  in  the  sample. 
g/,(z)  is  taken  as  a  phenomenological  parameter,  in  part 
because  the  degree  of  alloying  between  the  nominally  pure 
InAs  dots  and  GaAs  barrier  is  unknown  [19]. 

If  we  assume  that  the  g  factors  for  the  holes  localized 
in  each  dot  are  the  same,  we  get  gB(F)  =  gh  +  2 abgl2, 
where  g12  =  (l|g/,(z)|2)  gives  the  contribution  from  the 
amplitude  of  the  wave  function  in  the  barrier.  In  the  case  of 
the  X$  singlet  [shown  in  Figs.  3(c)  and  4(c)],  the  initial 
states  are  split  by  ge,  so  the  total  transition  g  factor  is  given 
by 

gBAF)=ge  +  gBh(F)  =  gT+  2thgn  . ,  (3) 

V e2d\F-Fh)2  +  4tl 

where  the  explicit  form  for  2 ab  determines  the  line  shape 
centered  around  the  anticrossing  point  Fh. 

The  lower  black  line  in  Fig.  5  is  obtained  by  fitting 
Eq.  (3)  to  XX Using  the  measured  values  of  th 
(0.86  meV)  and  Fh  (63.3  kV/cm)  we  find  g12  =  1.65. 
The  barrier  contribution  is  positive,  like  the  heavy  hole  g 
factor  in  bulk  GaAs  (—1.05)  [20,21].  Because  gT  and  g12 
have  opposite  sign,  the  splitting  reaches  a  minimum  at  Fh, 
where  the  amplitude  in  the  barrier  is  at  a  maximum.  The  g 
factor  at  the  minimum  is  —0.44. 


Electric  Field  (kV/cm) 

FIG.  5  (color  online).  Zeeman  splitting  and  corresponding  g 
factor,  gT,  as  a  function  of  electric  field  at  B  =  6  T. 
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Because  the  ±5/2  triplet  states  recombine  to  the  same 
bonding  orbital  of  the  final  hole  (Fig.  3),  the  model  predicts 
the  splitting  of  these  lines  as  a  function  of  electric  field  with 
no  additional  fitting.  As  shown  in  Fig.  4c,  the  ±5/2  triplet 
states  have  an  initial  state  splitting  of  ge  +  2gh.  The  tran¬ 
sition  g  factor  is  therefore  given  by  gf(F)  =  (ge  +  2 gh)  — 
{gi,  +  2 abgl2)  =  gT  ~  2 abgl2.  This  is  shown  by  the  upper 
black  line,  which  matches  AY/5^,  the  observed  splitting 
of  the  ±5/2  triplets.  The  maximum  splitting  corresponds 
to  a  g  factor  of  —4.23. 

To  apply  the  g-factor  model  to  the  bonding  orbital  of  the 
neutral  exciton  [lower  branch  Fig.  2(a)]  we  use  the  mea¬ 
sured  tunneling  energy  (tx o  =  0.58  meV)  and  anticrossing 
field  (Fx o  =  82.1  kV/cm).  The  lower  red  line  in  Fig.  5 
shows  the  fit  to  AYg,  the  Zeeman  splitting  of  PL  lines  from 
the  bonding  orbital.  We  find  g'l2  =  1-32.  The  electron-hole 
Coulomb  interaction  is  responsible  for  the  difference  in 
tunneling  energy  and  anticrossing  field  from  the  bare  hole 
case  and  is  also  likely  responsible  for  the  difference  be¬ 
tween  gl2  and  g j2.  The  g  factor  at  the  minimum  is  —0.59. 

Using  this  value  of  g\2.  the  model  immediately  explains 
the  increase  in  splitting  for  the  antibonding  orbital  [upper 
branch  Fig.  2(a)],  which  has  a  reduced  wave  function 
amplitude  within  the  barrier.  The  g  factor  for  the  antibond¬ 
ing  orbital  is  given  by  gj{F)  =  gT  —  2abg[2,  which  in¬ 
creases  in  magnitude  at  the  resonant  field  because  gT  is 
negative  and  g[2  is  positive.  This  is  shown  by  the  upper  red 
line  in  Fig.  5,  which  matches  AY®,  the  Zeeman  splitting  for 
excitonic  recombination  from  the  antibonding  orbital.  The 
splitting  increases  to  a  maximum  (g  factor  —3.35)  at  the 
anticrossing  point.  The  antibonding  transitions  for  Y+ 
show  similar  behavior,  but  are  too  weak  to  obtain  full 
resonance  curves. 

The  overall  agreement  between  the  model  and  experi¬ 
mental  data  is  quite  good.  There  are  some  minor  discrep¬ 
ancies,  which  highlight  the  need  for  a  detailed  theory, 
possibly  requiring  inclusion  of  excited  states  [22]. 
However,  the  agreement  of  the  data  with  the  resonance 
linewidths  calculated  using  independently  measured  val¬ 
ues  of  tx»  and  th  is  strong  confirmation  that  the  # -factor 
dependence  does  arise  from  the  formation  of  bonding  and 
antibonding  orbitals.  For  X+  the  g-factor  resonance  arises 
from  the  wavefunction  of  the  single  hole  in  the  final  state, 
while  for  Y°  the  orbital  wavefunction  of  the  hole  is  influ¬ 
enced  by  the  additional  electron. 

We  also  studied  samples  in  which  electrons  tunnel 
through  the  barrier,  with  an  anticrossing  energy 
(—2.3  meV)  comparable  to  that  of  the  hole-tunneling  sam¬ 
ple  presented  here  (—1.7  meV).  This  requires  a  thicker 
barrier  (10  nm)  because  of  the  smaller  electron  effective 
mass.  The  electron  wave  function  amplitudes  in  the  barrier 
should  be  at  least  as  large  as  the  hole-tunneling  case. 
However,  the  electron  g  factor  in  bulk  GaAs  (—0.44) 
[23]  is  similar  to  the  electron  g  factor  in  InAs  quantum 
dots  (—0.6)  and  according  to  our  model  the  contribution 
from  the  barrier  should  not  significantly  change  the  elec¬ 


tron  g  factor.  We  see  no  electric  field  dependence  of  the  g 
factor  in  these  electron  anticrossing  samples.  By  adding 
aluminum  to  the  barrier,  the  resonant  contribution  to  the 
electron  g  factor  could  be  enhanced. 

We  have  presented  a  resonant  change  in  g  factor  as  a 
function  of  electric  field  for  the  molecular  spin  states  of 
QDMs  with  a  thin  tunnel  barrier.  By  studying  single 
QDMs,  we  are  able  to  identify  the  individual  molecular 
spin  states  and  the  different  resonant  behavior  of  their  g 
factors.  The  results  suggest  that  design  of  molecular  spin 
states  and  tunnel  resonances  may  provide  new  opportuni¬ 
ties  for  combining  optical  and  electrical  control  of  confined 
spins. 
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